Upward shower rates at neutrino telescopes directly determine the neutrino flux 
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We show that the rate for upward showers from an isotropic cosmic neutrino flux at neutrino 
telescopes like IceCube is independent of the neutrino-nucleon cross section. For bins that span 
a relatively narrow range in energy, neither scaling the cross section, nor changing its power-law 
energy behavior affects the upward shower rate, which depends only on the flux. The neutrino 
flux can be completely known since its spectral shape can be determined by comparing the rates in 
neighboring bins. We also show that the downward shower rate varies linearly with cross section with 
a proportionality constant determined by the energy-dependence of the cross section, independent 
of the power-law behavior of the flux. The normalization and energy dependence of the cross section 
can be found by comparing the downward rates in neighboring bins. 



Introduction. Neutrino telescopes looking for high 
energy cosmic neutrinos are already ruling out jl| more 
aggressive flux models from diffuse sources Q and of 
cosmogenic origin 3]. The limits placed are within the 
framework of the Standard Model (SM). However, it is 
conceivable that new physics manifests itself in neutrino- 
nucleon interactions at high energies, enhancing the in- 
terest in and difficulty of studying high energy neutri- 
nos. The task of simultaneously determining the cosmic 
flux and the high energy cross section is a frustrating 
task since laboratory sources to examine neutrino inter- 
actions in controlled experiments at the needed energies 
are out of the question. Astrophysical fluxes cannot be 
tweaked in the control room, and the highest energies 
are beyond the reach of even the Large Hadron Collider. 
The challenge to find observables that lead to efficient 
extraction of the ffuxes and interactions is the motiva- 
tion for this paper, as it has been for a number already 
in print ^, 5, 6, 7]. 

We investigate the effect on shower rates of varying the 
power-law energy dependence of the cosmic neutrino flux 
and of the cross section in addition to the overall strength 
of the cross section. We illustrate how to extract com- 
plete knowledge (normalization and energy dependence) 
of the cosmic flux and the neutrino-nucleon cross section. 

Methodology. Using reasoning similar to that of 
Ref. 6], the rate for showers in a volume detector cen- 
tered at depth d can be written as 



dn (R - d)Xd 



\e-i{ei -1){^ + I]dl, (1) 



where ^ is the differential flux, Ap is the area projected 
against the neutrino direction, R is the radius of the 
earth, s is the scale of the detector size, I is the chord 
length traversed by a neutrino, and = \/ d(2R — d) 
is the horizontal distance from the detector to the "hori- 
zon" . Since the nucleon density n varies along a chord, 
we use an attentuation length A = l/{n(7t) that is aver- 
aged over I; at is the total interaction cross section. To 
emphasize that the interaction length at the detector is 
not averaged over a chord segment, we denote it as Xd- 
For up (down) going events, the integration limits are 



a = Lh, b = 2R — d {a = d, b = Lh)- The nadir angle 
that divides upward from downward events is therefore 
90°, as measured from the detector at depth d. We find 
that changing the definition of upward versus downward 
events by ±3° has no effect on our conclusion that the 
upward event rate is determined by the flux alone. 

While Eq. ^ accounts for the depth of the detector 
and the Earth's density profile [8.], it does not include 
downscattering or neutrino regeneration effects, which 
however, are typically of order 10% and irrelevant for 
our considerations. The full propagation requires the so- 
lution to coupled integro-differential equations, and, as 
shown in Ref. Q , the numerical calculations support the 
conclusions one draws from the analytic result. 

We employ the description of physics models intro- 
duced in Ref. @. The strength of charged current (c) 
and neutral current (n) interactions are parameterized 
by oic = ac/<J^^ and q;„ — (Tn/crf^, respectively, and are 
assumed to be energy-independent within a bin. The SM 
cross sections correspond to (pic, an) — {fcfn), where 
Ti = a^^/a^^. We consider new physics that scales 
the SM charged current and neutral current cross sec- 
tions by the same amount, (q;c,q;„) = (arc,Q;r„), with 
inelasticity the same as in the SM. We allow a to take 
values from 0.2 to 10. We also allow for changes in 
the power law behavior of the SM, which is approxi- 
mately (TsM{Ey) E^-'^^ [§]. We adopt a simple form 
(T(£:^j - £;0.36+/3 ^j^j^ p ranging from -0.36 to 1, and 
normalize so that /3 = reproduces the rate for a = 1} 

We consider neutrino fluxes that fall with 
ergy as E'^'^ , with 7 = 1,2,3 and normalize 
event rates at a = 1 to that for an isotropic 
tal flux of neutrinos and antineutrinos 4^ = 

. , dil 

10'8(£:^/GeV)-2 (cm^.s.sr.GeV)-! ^ with an assumed 
flavor ratio at Earth of 1:1:1; the Waxman-Bahcall 
(WB) flux. We integrate rates over bins of width 
Alogio(S,/GeV) = 0.5. 
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^ Note that power law growths of E^~^-^ are typical of inelas- 
tic, hadronic shower dominated cross sections from large extra 
dimension physics in the energy bins under study lid . 
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FIG. 1: Event rates (in three representative energy bins) 
vs. a, for fluxes falling with energy as E~~' . The plots are 
normalized at q = 1 to the WB event rate. 
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FIG. 2: Event rates vs. a, for fluxes falling with energy as 
E^'' , with their true normalization for the middle energy bin, 
10'' <E^< 10^-^ GeV. 



Dependence of shower rates on a and 7. We 

first set f3 — 0. In Fig. [H we plot the rate for up- 
ward showers in energy bins spaced equally in logarithm, 
lO^-s < < 10^ GeV, 10^ < < lO^-^ GeV and 
lO'^-S < E^ < 10^ GeV, as a function of a for 7 = 1,2 
and 3. For convenience in comparing the shape of the 
functional dependence of rates on a for different input 
fiuxes, we have normalized the 7=1 and 3 curves to 
the WB value at a = 1, in each energy bin. The result- 
ing shapes are nearly identical, independent of 7. For 
all three fluxes, the variation of the rate is a factor 2 or 
less over the factor 50 change in cross section strength. 
The variation is less than 50% for a between 0.5 and 5. 
The flux determined by using the SM cross section in the 
analysis is the correct one, since any other cross section 
with P = would give the same value. 

In Fig. [2] we show the curves for the middle bin with 
their true normalization, since the curves in each energy 
bin in Fig. [1] are normalized to the value appropriate to 
the WB flux with the SM cross section for convenience 
in comparing the shapes. Putting the information from 
the two figures together, we see clearly that the spectral 
index of the flux affects the rate in a bin, but it does not 
affect the independence of the rate from a. 

Dependence of shower rates on a, (3 and 7. We 
now explore the effect of varying a, [3 and 7 simul- 
taneously. In Fig. [3] we show how the upward event 
rates in our three representative bins are affected for 
/9 = —0.36,0,0.5 and 1.0. We only present the results 
for the WB fiux (7 — 2) because for a given /3, choosing 
7 = 1,3 makes a barely perceptible difference. Allowing 
(3 to range from —0.36 to 1 produces a shift in shape of 
the curves, but the values stay within the same 50% range 
around the SM value for any value of a. Neither energy 
dependence nor normalization of the cross section appre- 
ciably modifies the up shower event rate in our "fiducial" 
bins. Changing the energy dependence of the flux also 
makes little difference. The upward shower rate per bin is 
determined by the flux alone. The energy dependence of 
the flux can be determined by comparing rates in neigh- 
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FIG. 3: Upward shower event rates vs. a for the WB flux for 
four values of /3. The curves are normalized such that the SM 
rate is reproduced for /3 = and a = 1. 



boring energy bins, as we demonstrate later. 

The down shower rates as a function of a for different (3 
and 7 are shown in Fig. |4l The expected proportionality 
to a 6] is clear, with values of the rates determined by 
/3. We see that the event rate for /3 = is unaffected by 
changing 7, since the dashed line, the open circles and 
the plus signs lie on top of each other. 

Flux and cross section from data. Given the en- 
ergy dependence of the flux, the measured rate of upward 
showers in an energy bin will determine the normaliza- 
tion of the flux. It is the only other quantity, apart from 
the known acceptance and earth density proflle, that the 
rate in a bin depends upon over a wide range of model 
assumptions. To obtain the energy dependence of the 
flux, one compares the rates in neighboring bins to ex- 
tract 7 in that range of energy. To the extent that the 
neighboring bins which span the energy ranges Ei-i — Ei 
and Ei — £"^+1, are approximated by the same power law 
in energy, and the acceptance, is roughly constant in 
each bin, the ratio of upward shower rates in neighboring 
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FIG. 4: Similar to Fig. |3] but for downward showers only in 
the energy bin 10^ < < lO''-^ GeV. The squares almost 
on top of the dot-dashed line is the SM rate (/3 = 0) rescaled 
to agree with the /3 = 1 rate at a — 0.2. The slight slump of 
the /3 = 1 curve relative to the SM indicates that absorption 
at the highest cross section values is starting to appear. Note 
that the event rate for (3 — is unaffected by changing 7. 



bin[£,_i, E^]/hm[E^, E^+l] 
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TABLE I: Ratio of upward shower rates in adjacent bins for 
the WB flux from a simulation of the IceCube detector. 

bins is for 7^1, 

r{i,i + 1) ^7^-1 _ E]-^ \E^-lJ A^^^+i ' 

(2) 

where the rate in a bin is denoted by T(i,i + 1), 
the acceptance by and the average density ratio 

PEarth/pice Hiet by the incoming flux contributing to 
events in a bin by pi^i+i. For the special case 7=1, 
the explicitly energy dependent factors are replaced by 
ln{Ei / Ei-i) / ln{Ei+i / Ei) (When the bins are equally 
spaced, the ratios will all be 1.) In the table, we show 
the result of this simple analysis for simulated up shower 
data with the IceCube effective volume and the 
WB flux. We list the ratio of rates in neighboring bins 
from simulated data of upward showers in an "IceCube - 
like" volume detector. The actual values are taken for the 
SM cross section, but changing the cross section changes 
the ratios very little, as we have shown. Equation 2 im- 
plies that for all bins, the values for the 7 = 1,2 and 
3 flux spectra are 10°, 10° '^ and 10^, respectively, from 
the explicit energy factor, with reductions of order 10% 
coming from the acceptance and density factors. Clearly 
the values are consistent only for 7 = 2, which is the 
spectrum used to generate the data. 

A graphical representation is shown in Fig. [5l We use 
the up shower event rates generated by the SM cross 
section and WB flux in our representative bins to extract 
the overall coefficients, F-^, of the F^E~'^ flux spectra for 
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FIG. 5: Flux coefficient, F^, for fluxes falling with energy 
as EiJ~' . To display on the same plot, we multiplied the co- 
efficients by 10^^, 10*, and 1. The event rate data used in 
the analysis is generated from the WB flux model, which has 
7 = 2. Five energy bins from lO'''^ to 10^ GeV of width 
Alog^Q(_Ey/GeV) = 0.5 are chosen for illustration. The co- 
efficient must be the same in each bin to be consistent with 
data. Clearly, the WB case is reproduced while the 7 = 1 and 
3 cases are strongly excluded. 

7 = 1, 2 and 3. Factors appropriate to their different 
dimensions were applied to display them all on the same 
graph. With these simple input models, the data choose 
the one whose coefficient is the same in each bin. As 
shown in the figure, the 7=1 and 3 cases are easily 
excluded. As indicated by our earlier analysis, this result 
will be the same even for widely different cross sections. 

The fiux in a given bin can be estimated from Eq. p]). 
given the value of the effective volume T4/ f — sx Ap, Tup 
for the bin, the Earth's radius and the ratio of the average 
density of earth to density of target (ice). The average 
density the incoming neutrinos "see" depends mildly on 
energy. In the bin between 10^ GeV and 10^'^ GeV, sim- 
ulation gives ~ 0.3 events/year for cross sections with a 
between 0.2 and 10, with the normalization of Ref. pTl |. 
We use Veff ~ 2 km'^ [T2| and average density ratio 
PEarth/ Pice ~ 2 — 3 for the crust/mantel segment that 
dominates the upcoming flux, and flnd the flux normal- 
ization to be 6 — 9 • 10^^^ (cm^- s • sr)~^, which agrees 
roughly with the input normalization. 

Extracting the differential flux from data. Ex- 
panding on this question of flux measurement, we next 
estimate the differential flux index directly from data by 
the following prescription. Given data for upward shower 
rates in several adjoining energy bins, we refer again to 
Eq. Il]) to motivate the deflnition, for 7 7^ 1, 

d<j)iE,) _ R r(i - 1, vJ,X^-\/2) - » + 1) vJl'(CAli) 

(3) 

where the j ± 1/2 notation indicates that the quantity 
is evaluated at the midpoint (in logarithm) of the rele- 
vant bin. The energy assigned to the flux is between two 
adjoining bins for plotting purposes. The flux index is 
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FIG. 6: The squares show the output values of the differential 
flux extracted by the procedure of Eq. (|3]) (except for 7 ~ f ), 
from the numbers of up shower events calculated in our model. 
The event rate data that are used in the analysis are generated 
from flux models with 7 = 1, 2 and 3, all normalized to the 
WB flux at f 0^ GeV, and the cross sections have f3 = and 
a = 1, 1 and 6, respectively. Five energy bins from 10^'^ to 
iO^ GeV of width A\og^g{E^/GeV) = 0.5 are chosen for the 
analysis, which then produces flux values at four energies. 

found by fitting the points obtained from Eq. ([3]). The 
normahzation can then be obtained from the event rate 
in a particular bin. If the observed rates are nearly the 
same in each bin, 7 ~ 1 applies and Eq. ([3]) is irrelevant. 
A fit near 7 = 1 can be found and the normalization 
determined from any bin in this range. Applying this 
prescription to our analytic model for three fluxes, we 
obtain the values for the differential fluxes shown by the 
squares in Fig. [S) The lines are the input fluxes nor- 
malized to the WB flux at 10^ GeV. The procedure of 
Eq. ((21) is obviously successful in extracting the original 
flux. We have checked that the procedure is effective for 
other flux and cross section combinations, as expected 



from the results shown in Figs.[Tl[2]and[3l 

The ratio of downward events to upward events in a 
given energy bin is potentially a clean way to extract 
information about the total cross section in that bin 
[Tsj . Using Eq. ([T]) in the usual energy range, with average 
density 2 — 3 • Na, and a ratio of down to up events ~ 7, 
we find cr = 3 — 5 ■ 10~^^ cm^ x a. This agrees nicely with 
the SM value at a = 1 Q in this energy interval. 

Summary. Our study confirms the linear dependence 
of downward events on the product of cross section and 
flux values and establishes the robustness of the linear 
dependence of upward shower event rates on flux alone. 
Comparing rates in neighboring bins yields the energy 
and normalization of flux and cross section. We have 
focused on the systematics of the analysis by using a 
model which captures the main features of the upcom- 
ing shower events. It has the essential elements of a full 
simulation, but of course does not pretend to provide the 
detail necessary for analysis of real data. We have not 
addressed the tough question of reconstructing neutrino 
energy, the variable we adopt directly here, from the en- 
ergy measured in a shower. 

Our message here is that distinguishing upcoming 
showers from downgoing ones, with reasonable energy 
resolution can have a big payoff in scientific discovery. If 
the charged current events are proportional to the total 
cross section, the discussion we presented applies equally 
well to them. As shown here, the upcoming showers pro- 
vide a remarkably clean way to isolate the flux indepen- 
dently of cross section, over a wide range of energies. 
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